Measuring errors in single qubit rotations by pulsed electron paramagnetic resonance 
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The ability to measure and reduce systematic errors in single-qubit logic gates is crucial when eval- 
uating quantum computing implementations. We describe pulsed electron paramagnetic resonance 
(EPR) sequences that can be used to measure precisely even small systematic errors in rotations 
of electron-spin-based qubits. Using these sequences we obtain values for errors in rotation angle 
and axis for single-qubit rotations using a commercial EPR spectrometer. We conclude that errors 
in qubit operations by pulsed EPR are not limiting factors in the implementation of electron-spin 
based quantum computers. 
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Pulsed magnetic resonance methods have provided a 
useful playground in which to test different aspects of 
quantum computation. Single-qubit operations (spin ro- 
tations) can be conveniently performed using classical ra- 
dio frequency (RF) pulses whilst two-qubit operations 
can be naturally realized through exchange or dipolar in- 
teractions. Exploiting these advantages, the largest-scale 
quantum computations to date have been demonstrated 
using nuclear magnetic resonance (NMR) in liquid solu- 
tion [1, 2]. Electron paramagnetic resonance (EPR) offers 
many parallels with NMR, along with the key advantage 
that the electron gyromagnetic ratio is of the order of a 
thousand times larger. Thus pure ground states are ex- 
perimentally accessible in EPR, avoiding the scalability 
issues surrounding NMR implementations [3] . 

For this reason, EPR has become a key element in sev- 
eral solid state quantum information processing (QIP) 
proposals [4-8]. The practicality of these proposals is 
critically dependent on the errors that are inherent in 
pulsed magnetic resonance experiments. In particular, 
while the decoherence time (T2) is generally quoted as 
the ultimate figure of merit for qubit implementations, 
understanding and minimizing the systematic errors in- 
herent in qubit manipulations is equally important. A 
number of general approaches to tackling different classes 
of systematic error in qubit rotations, employing com- 
posite rotation sequences, have been proposed [9, 10]. 
However, before these approaches can be exploited prac- 
tically, it is necessary to characterize (and reduce) the 
errors that are associated with a single rotation pulse. 

There are two principal types of systematic error asso- 
ciated with a single-qubit rotation: rotation angle error 
and rotation axis error. In magnetic resonance experi- 
ments, rotation angle errors arise from an uncertainty in 
the Rabi oscillation period, associated with uncertainty 
in the magnitude and duration of the applied RF pulse. 
On the other hand, rotation axis errors arise from un- 
certainty in direction of the RF magnetic field in the 



transverse plane of the rotating frame. 

In this Letter, we show how multi-pulse sequences can 
be applied to measure precisely these two classes of er- 
ror. These sequences include some originally developed 
for NMR, as well as a novel sequence developed specif- 
ically for quantifying phase errors. In each case the se- 
quence amplifies the errors, so that even small errors are 
detectable. We find that when comparing a commercial 
EPR spectrometer to an NMR system which has been 
optimised for high-fidelity qubit operations, rotation an- 
gle errors are as good, whilst phase errors in EPR are 
worse by as much as an order of magnitude. 

The paramagnetic species used here to perform er- 
ror measurements in pulsed EPR is i-NCeo (also known 
as N@Cqo), consisting of an isolated nitrogen atom in 
the 4 S 3 / 2 electronic state incarcerated by a Cgo fullerene 
cage. It is an ideal system for these measurements be- 
cause of its extremely narrow EPR linewidth and long 
relaxation time in liquid solution [11, 12]. T2 has been 
measured to be 80 (is at room temperature, rising to 
240 fis at 170 K [13]. 

The production and subsequent purification of i-NCgo 
is described elsewhere [14]. High-purity i-NCgo powder 
was dissolved in CS2 to a final concentration of 10 15 /cm 3 , 
freeze-pumped in three cycles to remove oxygen, and fi- 
nally sealed in a quartz EPR tube. Samples were 0.7- 
1.4 cm long, and contained approximately 5 • 10 13 i- 
NCqo spins. Pulsed EPR measurements were done at 
190 K using an X-band Bruker Elexsys580e spectrometer, 
equipped with a nitrogen- flow cryostat ( Janis Research) . 

i-NCeo has electron spin S = 3/2 coupled to the 14 N 
nuclear spin 1 = 1. The EPR spectrum consists of three 
lines centered at electron g-factor g = 2.003 and split by 
14 N hyperfine interaction a = 0.56 mT in CS2 [15]. How- 
ever, all pulsed EPR experiments discussed below were 
done using the center line in the EPR triplet, correspond- 
ing to the 14 N nuclear spin projection Mi = 0. Given the 
small isotropic hyperfine coupling in i-NCgo, the transi- 
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FIG. 1: Rabi oscillations for i-NCeo in CS 2 at 190 K. 



tions with simultaneous flip of electron and nuclear spins 
are largely forbidden and therefore the evolution of elec- 
tron spin can be treated individually for each nuclear spin 
manifold. The Bloch sphere is a useful aid to visualise 
the action of the pulse sequences described below [16]. 
While this model maps conveniently the magnetization 
evolution of a single S spin system, the evolution of a 
coupled S = 3/2, 1 = 1 system is harder to visualize. It 
can be shown that the evolution of the electron spin in 
the Mi = manifold is largely unaffected by the presence 
of the hypcrfinc coupling and is adequately described by 
the classical vector model [13]. 

In an EPR experiment, a qubit rotation is achieved 
by applying an on-resonancc microwave pulse of con- 
trolled power and duration. With a pulsed magnetic field 
strength B\ and a pulse duration t the rotation angle is: 



9 = g^sBit/h. 



(1) 



Possible off-resonance effects are neglected here because 
the EPR linewidth for i-NCeo is much smaller than the 
excitation width of the RF pulses used in our experi- 
ments. According to Eqn. 1, rotation angle errors may 
arise from either pulse duration errors (which can be as- 
sumed uniform throughout the sample), or errors in the 
magnitude of B\ (which can vary across the sample de- 
pending on the homogeneity of the EPR cavity mode). 
The hardware of the EPR spectrometer limits pulse du- 
ration to a 2 ns resolution and pulse power to 0.1 dB 
resolution. These together can contribute of the order of 
0.5% to the rotation angle error for a typical 32 ns pulse. 

In EPR, measurements are made of the magnetisation 
of the ensemble in the x — y plane, for a static magnetic 
field applied along z. Hence, a simple measurement of 
the rotation angle error could be made after a nominal 
7r/2 rotation, which in practice is it/2 + 5 where S is 
the error. The measured signal is proportional to cos(5). 
However, this approach is unsatisfactory for two reasons: 
for a good measurement, a reference (an ideal rotation) 
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FIG. 2: Comparison of the echo signal decays in (a) the CP 
and (b) the CPMG pulse sequences. The narrow spikes cor- 
respond to the applied tt pulses and the echo signals are seen 
in between, (c) Decay of echo magnitudes for each sequence 
(CP: empty, CPMG: filled). The solid grey line is a fit to a 
simple exponential yielding T2 = 190 /is; the dashed grey line 
is a fit yielding rotation angle error parameters [17]. 



is needed with which to compare the imperfect rotation; 
and the cos((5) term depends on 6 only to second order. 

Another method involves applying a long RF pulse 
and observing Rabi oscillations over a number of peri- 
ods. The respective experiment for z-NCeo is shown in 
Fig. 1, demonstrating Rabi oscillations whose amplitudes 
decay at long pulse durations (over 80 oscillations were 
seen) . This decay is caused by inhomogeneity of B\ fields 
in the EPR cavity (spins are rotated with slightly differ- 
ent Rabi frequencies and therefore gradually lose coher- 
ence), as well as effects such as B field inhomogeneity 
and the fact that the output power and phase from the 
microwave amplifier vary at long pulse lengths. In or- 
der to distill the rotation angle error from these other 
effects it is necessary to use more sophisticated methods 
for error measurement. 

A better approach uses sequences of many pulses to 
compound the error of a single pulse. A suitable se- 
quence, CP (Carr-Purcell [18]), consists of a 7r/2 pulse 
followed by a series of rcfocusing tt pulses, i.e. ir/2 x — 
(t — ir x — t)„ in which pulse rotation angle imperfec- 
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FIG. 3: Evolution of the magnetisation vector in the rotating 
frame during a SPAM experiment. The ideal y and x axes 
are associated with in-phase and quadrature channels of the 
quadrature detector, respectively. tt/2 x and tt x pulses are 
assumed to be along the x-axis; n y pulses deviate from the 
ideal y-axis by the error angle S and are oriented along y' . 
The magnetisation vector (bold) is shown (a) immediately 
after the ir/2 x pulse, and at times of echo formation after (b) 
the first n y pulse, (c) the first tt x pulse, (d) the second ir y 
pulse. Alternating the refocussing pulses between x and y' in 
SPAM results in an accumulation of error in the phase of the 
echo signal. 



tions are additive [19]. The typical exponential decay of 
the echo is therefore further attenuated by the cumula- 
tive rotation angle error, such that the echo amplitude 
decays with a time constant shorter than Ti. Meiboom 
and Gill proposed a modification of the CP sequence, 
termed CPMG [20], which compensates for pulse length 
errors by applying the refocussing tt pulses around the 
y— axis. 

Figs. 2(a) and (b) compare the trains of echo signals 
observed in the CP and CPMG experiments for i-NCgo- 
The decay of the CPMG echo magnitudes is fitted to 
a simple exponential to obtain the transverse relaxation 
time Ti. In determining the decay of echo magnitudes 
due to rotation angle errors in the CP sequence, we as- 
sume the flip angle error follows a Gaussian distribution 
with mean 5o (to account for pulse duration errors) and 
standard deviation as (to account for inhomogcncity in 
the oscillatory magnetic field strength) [17]. From the fit 
we obtain as = 18° in every 180° rotation, or approxi- 
mately 10%. This figure is consistent with the expected 
inhomogcncity in the applied B\ field, and was seen to 
vary by altering the dimensions of the sample. 

Having determined the rotation angle errors, we now 
turn to rotation axis errors. The Bruker spectrometer 
used in this work offers four independent pulse-forming 
channels, each supplied with uncalibrated analogue phase 
shifter. Our initial goal was to set microwave phases 
in two pulse-forming channels to be orthogonal to each 
other, which is to orient the B\ field in one channel along 
the x-axis and the other along the y-axis. The ideal x 
and y-axes in the rotating frame are defined with respect 
to the phase of quadrature detection channels, therefore, 
before perfecting the phase setting of the pulse-forming 
channels we first examined the orthogonality between the 
two channels of the quadrature detector. This was done 
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FIG. 4: (a) The echo train from SPAM showing a leakage of 
signal between the y (in-phase: top trace) and x (quadrature: 
bottom trace) detection channels resulting from phase error 
between nominal n x and ir y pulses (traces are offset for clar- 
ity), (b) Magnitudes of successive echoes from (a) fitted to 
functions described in the text yield 6 = (10.3 ± 0.5)°. 



by applying a slightly off- resonance 7r/2 pulse and ob- 
serving the resulting FID oscillations from each detection 
channel. A numerical fit revealed an angle of 89.3±1° be- 
tween the nominally orthogonal "real" and "imaginary" 
detection channels. 

Traditionally, the phase of a pulse-forming channel is 
adjusted by applying a simple 7r/2 rotation about each 
channel and observing the FID signal. The inherent im- 
precision of this approach is not a serious problem in 
traditional EPR applications in which only a few pulses 
are applied, but is potentially devastating for the fidelity 
of a qubit state in a multiple-pulse computation. 

In order to measure the orthogonality of the rotation 
axes in two pulse-forming channels with a high preci- 
sion we designed a pulse sequence to accumulate phase 
errors, by analogy with the way that the CP sequence 
above accumulates flip angle error. This sequence, called 
Sequence for Phase-error AMplification (SPAM) is given 
in Eqn. 2 and illustrated in Fig. 3. 



ir/2 x (t TTy T T 71 x 7~) r 



(2) 



After n cycles of this sequence with a phase error <5 (i.e. 
the phase of the 7r y pulses is actually n/2 + S with re- 
spect to the tt x pulses) the echo vector points along 
(cos(2n5), sin(2n<5), 0). Thus, a non-zero 8 in this sc- 
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quence results in an accumulating "leakage" of the echo 
amplitude into the orthogonal detection channel, with 
only a second-order sensitivity to the pulse length errors 
described above. 

Fig. 4(a) shows the echo train obtained in a SPAM se- 
quence for Z-NC60- For this experiment a phase error of 
approximately 10° between the tt x and tt v pulses was in- 
tentionally introduced by observing the FID signals. The 
measured echo magnitudes for the real and imaginary 
channels are described by cos (2n5) and sin (2nS) respec- 
tively, with a exponential decay envelope exp (— fair jT-i). 
n is the number of SPAM cycles (each comprising two 
refocusing pulses) and Ti is taken from a CPMG exper- 
iment described above. Fig. 4(b) shows the echo ampli- 
tudes extracted from (a). The fit yields 5 = (10.3 ±0.5)°, 
which is close to the intended phase error. This demon- 
strates that SPAM is a useful way of both measuring the 
phase error in nominally orthogonal channels and accu- 
rately setting arbitrary phases between channels. 

Orthogonality of the ir x and ir y channels was then op- 
timised by traditional procedures (by observing the FID 
as described above), and a phase error of S = (1.5 ±0.3)° 
was measured from a SPAM sequence. This provides a 
measure of the phase error typical in conventional pulsed 
EPR experiments. The fact that it can be measured to 
a precision of about 0.3° using a SPAM sequence allows 
us to reduce it substantially further. Careful optimisa- 
tion of the phase error using the SPAM sequence yielded 
6 = (0.3 ±0.1)°. 

In summary, we find that although commercial pulsed 
EPR spectrometers have not been designed with high 
precision operations and multi-pulse sequences in mind, 
the Brukcr machine performs very well. The rotation 
angle error of about 10% is almost entirely due to inho- 
mogeneity of the oscillatory magnetic field, and is com- 
parable with those typically observed in NMR [10]. The 
rotation axis (phase) error of about 0.3° is also compa- 
rable with those typically encountered in NMR (though 
worse than in optimised quantum computing NMR spec- 
trometers). 

The most significant strength of the SPAM methodol- 
ogy is that it provides a method of setting relative phases 
between channels with very high precision. This makes it 
possible to exploit techniques developed within the con- 
text of NMR for applying sequences of pulses at various 
phases that correct for rotation angle errors [10, 21, 22]. 
Such a sequence can straightforwardly reduce a rotation 
angle error of order e to e 6 . Using this approach, our 
results imply that rotation angle errors can be reduced 
the order of 10~ 6 , well within the threshold of 10~ 4 often 
cited for fault-tolerant quantum computation [23]. This 
is the subject of a subsequent paper. 

We have shown that the technology exists to transfer 
the quantum information processing methodology that 
has been developed in the context of NMR to EPR, 
thereby overcoming the scaling limitations associated 



with NMR. We conclude that errors in qubit operations 
do not restrict the viability of an EPR-based quantum 
computer. Furthermore, we demonstrate the successful 
application of long pulse sequences (necessary for run- 
ning quantum algorithms) to an EPR qubit candidate, 
the z-NC6o molecule. Finally, we have demonstrated a 
set of pulse sequences that can be used to amplify and 
measure precisely the effect of rotation angle and axis 
errors in any NMR or EPR pulsed magnetic resonance 
spectrometer. 
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SUPPORTING MATERIAL 

In the Carr-Purcell (CP) sequence, n/2 x — (r— ir x — r) n , 
the effect of rotation angle errors is complicated by the 
dispersion of spins in the x — y plane due to Bq field 
inhomogeneity. For example, those spins pointing along 
y when the refocussing ir x pulse is applied pick up the 
most error, whilst those pointing along x are unaffected. 
We assume a uniform distribution of phase, i.e. t >> T 2 *, 
and a Gaussian distribution of rotation angle error with 
mean So and standard deviation a. The echo magnitudes 
obey the following equation, after the nth pulse in the 
sequence: 



A C p{n) = f a m + y~] frfc exp ( cos (kS ) 



m—l \ k—1 



(1) 



2m Cm "+^-i C ' 2m _ 1 i/2 Cm n (2m-l) 

2m 1 > 



, 2m C m - fe "+ m - 1 C 2m _ 1 1 /2 Cm n (2m - 1) 

m 

(3) 

For no < 1, and assuming <5o = 0, this can be approxi- 
mated to Eqn. (4) below. 



A C p(n)=exp(——) (4) 



To fit the CP decay we use Acp(n) exp (— t/T 2 ), with 
T 2 taken from the CPMG fit. 



